Continuous vortex pumping into a spinor condensate with magnetic fields 
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We study the mechanisms and the Hmits of pumping vorticity into a spinor condensate through 
manipulations of magnetic (B-) fields. We discover a fundamental connection between the geomet- 
rical properties of the magnetic fields and the quantized circulation of magnetically trapped atoms, 
a result which generalizes several recent experimental and theoretical studies. The optimal proce- 
dures are devised that are capable of continuously increasing or decreasing a condensate's vorticity 
by repeating certain two step B-field manipulation protocols. We carry out detailed numerical sim- 
ulations that support the claim that our protocols are highly efficient, stable, and robust against 
small imperfections of all types. Our protocols can be implemented experimentally within current 
technologies. 
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I. INTRODUCTION 

A quantized vortex represents a hallmark of superflu- 
idity [l|. The vorticity in a superfluid clearly reveals 
the topological nature associated with the phase of a 
condensate's wave function and has been studied exten- 
sively both theoretically and experimentally 0, |2|] since 
the first success of atomic Bose-Einstein condensation. In 
the strongly rotating limit when the number of vortices 
significantly exceeds the number of atoms, researchers 
have focused on the possibilities of observing interes ting 
strongly correlated states 0, i, i, i, 0, 0, S llO, [HI, [H 
[H, [ij,ll5|,|16|]. 

Vortex states have been created experimentally relying 
on a variety of approaches: a spatially selective rotat- 
ing Rabi coupling to an auxiliary internal state [t?, ^S] , 
stirring inside [19] or at the edge [20] of a condensate 
with a focused laser beam, rotating a deformed trap |2l| , 
or simply merging several condensate pieces ^] into a 
single trap. Other approaches include direct phase engi- 
neering through the imprinting of a nontrivial topologi- 
cal phase to the wave function of a condensate [23.] us- 
ing an azimuthal optical phase plate or through ma- 
nipulating external magnetic (B-) fields [25l, [26l, [27l, [28j . 
Alternately, a condensate can gain angular momentum 
through its coupling to a Laguerre- Gaussian light beam 
pol , [so, [31', or other forms of gauge potentials [33] . 

The idea of phase imprinting from the nontrivial geo- 
metrical properties of external B-fields is especially illu- 
minating as it involves no forced rotation or essentially 
no spatial motion [25[. It has since been demonstrated 
in beautiful recent experiments [H, [35|, [s^. The orig- 
inal protocol [25I] involves an axially symmetric spinor 
condensate in the familiar B-field of an loffe-Pritchard 
trap (IPT) [37] consisting of a two dimensional (2D) 
quadrupole B-field augmented by a bias B-field along the 
symmetric z-axis. Upon the flipping of the bias field, the 
adiabatically following weak field seeking state then de- 
velops a nonzero vorticity in its spatial wavefunction [25]. 
Our aim is to develop this protocol into a pump mecha- 
nism capable of continuously changing the vorticity of a 



condensate. The naive approach of repeated flips of the 
z-bias B-field does not work because the second flip sim- 
ply undoes the vorticity gained in the first flip, returning 
the condensate to the initial state. 

Recently, Mottonen, et al [38], put forward an inter- 
esting idea that sought to break the time reversal sym- 
metry between the first and the second flips of the bias 
field. Their protocol starts with a 2D hexapole B-field in- 
stead of a quadrupole field for the first flip. The hexapole 
field is then turned off and replaced by a quadrupole 
one for the second flip. It is known, as in the origi- 
nal bias flip protocol [25|, that the first flip can gener- 
ate a proportionally higher vorticity if the quadrupole 
is replaced by a higher order 2D multipole such as a 
hexapole [s^. Repeating the two step protocols with 
hexapole and quadrupole B-fields in turns, Mottonen, et 
al.^ show that the vorticity increases by Ah per atom in 
the hexapole step and decreases hy 2h per atom during 
the quadrupole step [38] . Thus each cycle composed of a 
bias flip with a hexapole field followed by a second bias 
flip with a quadrupole field increases the net vorticity 
hy Aft — 2% = 2%. Higher vorticities are generated with 
repeated cycles. 

A related earlier discussion [3^ suggested an even sim- 
pler protocol of turning off the quadrupole field after the 
first flip, and then returning the axial bias field to its 
original direction. The return of the axial bias to the 
original direction is hoped to bring the system to the ini- 
tial configuration except for the vorticity gained from the 
first flip. It was previously noted that by turning on the 
quadrupole field and repeating the above protocol, a con- 
tinuous vortex pump is realized [39]. Unfortunately, no 
details were given on how the axial bias field is returned 
to the original direction [39] . Clearly it cannot be flipped 
back in the absence of the quadrupole field as the z-bias 
field remains aligned along the same axial direction. Sim- 
ply let it oscillate back with the axial B-field magnitude 
increasing from large negative to small negative, to zero, 
to small positive, and then finally to large positive val- 
ues does not constitute a flip. Both the projections of 
mechanical angular momentum Lz and the spin Fz{Mf) 
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are independently conserved quantities during the above 
process. 



II. OUR CONTINUOUS VORTEX PUMP 
PROTOCOL 

Our idea is physically intuitive and leads to a direct im- 
plementation for a continuous vortex pump with the in- 
troduction of an auxiliary transverse bias B-field. As in a 
Stern- Gerlach experiment, the Zeeman population distri- 
bution inside an experimental apparatus can be measured 
along any quantization direction of a strong reference B- 
field, provided the strong field is turned on adiabatically. 
When the transverse bias B-field is pulsed on, it provides 
a reference direction for the axial bias field to fiip back 
in the absence of the quadrupole field. We find it imper- 
ative to present this result because of the high interests 
stimulated by the recent work of Mottonen, et al ^] . 
Unlike the Ref. [s^ , our proposal requires only one set of 
current coils capable of generating one type of multipole 
B-field: be it a 2D quadrupole or a 2D hexapole. The 
time reversal symmetry is broken in our protocol with 
the extra transverse bias B-field whose presence allows 
for the 2;-bias field to rotate instead of simply oscillating 
back. 

As an illustration we consider a spin-1 condensate 
in an IPT, whose B-field is approximately B{x, z ) = 
B'{xx — yy) -\- BzZ near the origin. Isoshima et al. [25| 
first discussed pumping vorticity into a condensate with 
external B-fields. For a sufficiently large z-bias B-field, a 
condensate of F = 1 atoms adiabatically stays in the B- 
quantized \Mf = — 1, f)^ state. After an adiabatic fiip of 
the bias field, a vortical phase structure or a vortex state 
is imprinted into the condensate [26|, 127|, |28|, IM, l35|, l36| . 
The physics involved can be elucidated in terms of the 
symmetries and conservation laws of the model system 
[4^ or from the gauge potential [4l] due to the changing 
B-field. More specifically, it is the conservation of Dz 
(= Lz — Fz) that enables a straightforward understand- 
ing of the result of Isoshima, et al. [25], in the axially 
symmetric IPT [42|, |43| . Each single fiip of the bias field 
in the presence of a quadrupole field then imparts a 2Fh 
(for F = 1 here) phase winding [H, IH, [sB, [s^. Re- 
peated jyperations then constitute a continuous vortex 
pump [38]. As introduced above, our pump mechanism 
generalizes the original idea of Isoshima, et al. [25]. Fol- 
lowing the first flip, which imprints a 2h vortex to the 
adiabatic state, we turn on a transverse bias field B^x 
and then continue with a second bias fiip returning the 
system to the initial setup. A vortex pump then simply 
consists of repeated applications of the above manipula- 
tions to the three B-fields: the 2D quadrupole, the axial 
bias, and the transverse bias fields. 

The potential for high fidelity operation of our pump 
protocol is confirmed with numerical simulations for a 
spin-1 atomic condensate in external magnetic (B-) fields. 
An additional optical trap Vo = Mu'j_{x'^ + + A^z^)/2 



provides a permanent confinement during the B-field ma- 
nipulation. To avoid the energetic instability of the po- 
tential disintegration of a high vortex state into single 
vortices [50|, Ell, |52j] , we introduce an optical pinning plug 
Vp{p^z,(f)) = [/exp(— p^/2pq) that expels atoms away 
from the low B-field region to assure adiabaticity. The 
force of gravity is assumed to be opposite to the axial z- 
axis, which can be omitted in the approximate 2D treat- 
ment when the harmonic trap is taken to be pancake 
shaped with A 1. The z-dependence is frozen in the 
ground state of the axial harmonic oscillator (l){z). Thus 
the vorticity is coded in the phase structures of the two 
dimensional condensate wave function t). 

The condensate wave function thus be- 
comes ^(r, t) = t)^(z), where p = {x^y), 
(P{z) = (A/^ai)i/4e-^"'/2«i, and = ^fi/Muj^ 
is the length scale for the transverse harmonic trap. 
After integrating out the z coordinate, we obtain the 
effective 2D Gross-Pitaevskii equation 
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\MujI_{x^ + y'^) and V^^ = /7exp(-pV2p2) are the 2D 
optical trap and optical plug, respectively, n = 
is the 2D atom number density. The effective 2D spin- 
independent and spin-dependent interaction strengths 
are now characterized by Cq^^^ = 2{27rXy^'^h^{ao + 

2a2)/3Max, and 4^^^ = 2{27rXy/'^h^{a2 - ao)/3Ma±. 

In the local B-quantized representation, the Zeeman 
energy is diagonal, given by the Breit-Rabi formula 

Em, = ^ - Mpgif^iB - + Mf^ + e ■ (2) 

AEq is the hyperfine splitting, gj is the Lande factor for 
the nuclear spin /, pi is the nuclear magneton, and ^ is 
defined by ^ = {gipi + gjjiiB)B / AEq. Here gj is the 
Lande factor for the valence electron with total angular 
momentum J, and /i^ is the Bohr magneton. Our equa- 
tion ([1]) above requires Zeeman energy of an atom 
in the laboratory based z-axis quantization representa- 
tion that is nondiagonal and is given by = W ElA 
with U the unitary transformation connecting the z- and 
B-quantized representations. For the B-field of a IPT, 
the transformation matrix becomes the product of two 
rotations in spin space, i.e., U = e^^v^ e^^zfk ^ matrix 
form, 
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FIG. 1: (Color online), (a) The time dependence of the B- 
fields. In the first 30 ms, Bz decreases to the opposite of the 
initial value. The gradient increases from to the maximum 
in 15 ms, then decreases to zero. Next, the B-field is rotated 
back, while a uniform transverse bias in the x-axis direction 
is pulsed on for the latter 10 ms, keeping the + (B^)'^ 
a constant, (b) The time dependent fractional population 
riMp = Nmf/N in the Zeeman state. 



with the rotation angles and <p introduced through the 
parametrization of the field B = B[z cos ^+sin 0{x cos 
ysiiKp)]. Hf^^, is one of the nine elements of 
with spinor components ordered as Mp = 1,0,-1. 

Our numerical simulations assume a condensate of 
N = 10^, ^^Rb atoms in the F = 1 state and trapped as 
described above with A = 100. We take io± = 27r x 30 
Hz, U/h = 2 X 10^ Hz, and po = b pm. The B-field 
takes a form B = B' {xx — yy) -\- B^z -\- B^x with a typ- 
ical temporal evolution within one period as illustrated 
in Fig. [TJ During the first 15 ms, we decrease Bz and 
increase B^ according to the following easily programmed 
time dependence 



B\t) 



Bz(0)(e-8*/^i+2 - e-2)/(e2 - 
B'(Ti/2)(e2-e-8*/^i+2)/(e2 



(4) 

e-'), (5) 



with Ti = 30 ms. In the subsequent 15 ms, Bz is de- 
creased continuously as is the gradient B' . Both time 
dependence are relatively smooth and of an exponential 
type similar to the first 15 ms. The time dependence 
of the z-bias flip and the transverse x-bias pulse on are 
described by 



Bz(t) = Bz{0)cos[7t{1 
Blit) = 5^(0)sin[7r(l - 



-(t-T,)/T,)l (6) 
(t-Ti)/T2)], (7) 



with T2 = 10 ms. Due to the adiabatic nature of our 
proposed protocol, the quality of the final result or the 
fidelity of the intended vortex state do not sensitively 
depend on the details of the time dependence of the B- 
fields, provided they are reasonably smooth functions of 
time. In the shadow window of Fig. [TJ nothing really 
happens to the condensate spinor component population 



FIG. 2: (Color online). The temporal development of the 
condensate phase structures in a hexapole field Bh — B^^ [(x^ — 
y^)x — 2xyy] plus a transverse bias field BI^x. The upper row 
shows the Mp = 1 component at the end of the first flips, or 
in the middle of each cycle at 30, 70, 110, 150, and 190 ms, 
respectively; the lower row is for the Mf = — 1 state at 10 ms 
later with respect to the first row or at the end of each cycle 
at 40, 80, 120, 160, and 200 ms. The temporal dependence 
of the B-fields and other parameters are as in Fig. 1 of the 
supplementary material. The maximum for B^, Bz, and B^ 
are 1250 Gauss/cm^ 0.1, and 0.1 Gauss, respectively. White 
(black) color denotes phase — tt (tt). 



distribution as long as the B-fields change smoothly. Be- 
cause Bz is much larger than the field due to the gradient 
B^ within this window, the angle between the net B-field 
and the z-axis remains small even for a discontinuous 
change of the bias field or the quadruple gradient, caus- 
ing the atomic spin state to change smoothly. Depending 
on the net radial trap strength from the combined trap 
and the optical plug, the azimuthal vortical phase dis- 
tribution changes with the radius due to non-adiabatic 
effects with respect to the radial motional state. 

In all of our simulations, the initial state is obtained 
from the imaginary time propagation of the coupled 
Gross-Pitaevskii equations with an initial Gaussian wave 
function in the Mp = — 1 state setting C2 = and with no 
B-field as discussed in the supplementary material. The 
vortex pump protocol is then simulated in real time. For 
the case of a 2D quadrupole, the first flip of the bias Bz 
from 0.2 Gauss to —0.2 Gauss is carried out in the first 30 
ms, along with the increasing of the B-field gradient B^ to 
the maximum value of 10 Gauss/cm and followed with a 
gradual switch off. In the following 10 ms the bias field is 
rotated back around the 7/-axis in the x-z plane, aided by 
the transverse bias field along the x-axis. Similar steps 
are involved when the 2D quadrupole is replaced by a 2D 
hexapole. In Fig. [2]we show the temporal development of 
the phase structures for the Mf = 1 component after the 
first Bz flip and the Mp = — 1 component after the next 
flip back to the original direction for the first five repeti- 
tions. The total vorticity changes are impressive and con- 
firm conclusively our suggested protocol for a continuous 
vortex pump applied with a hexapole. For a quadrupole 
field, we find essentially the same quality operation, ex- 
cept that each repetition adds 2h units of vorticity as 
compared to Ah units for a hexapole (Fig. [2]). More gen- 
erally with a z-bias field along the symmetry axis of a 
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2/-th multipole 2D B-field B ex p^-^[cos{l(j))p - sin(/^)^], 
= Lz — {I — 1) X Fz is conserved for the symmet- 
ric ground condensate state [i^. A corresponding adia- 
batic flip of the bias field is thus capable of generating a 
2{l - l)Fh vortex. 

Before providing the general considerations that will 
lead to our discovery of the optimal vortex pump proto- 
cols based on manipulating external B-fields, we provide 
further numerical studies that demonstrate the efficiency, 
stability, and robustness of our vortex pump protocol. 

First, we discuss the relatively relaxed conditions on 
adiabaticity for the manipulation of the external B-field. 
Typically, the 2;-bias field is fiipped over during a time 
of 30 ms in our numerical simulations, although we find 
integer numbers of vorticity are still created provided the 
fiipping time is adjusted considerably (shorter or longer). 
Due to the symmetries and the corresponding conserva- 
tion laws on the dynamics of our model system; the con- 
servations of = Lz — Fz in an IPT, Lz + in a 3D 
quadrupole trap (QT) and Lz — {1 — l)Fz in a 21-th multi- 
pole 2D B-field, high fidelity operations are assured even 
with marginal conditions for adiabaticity. On careful ex- 
amination, we find that the aximuthal distributions of 
the phases are complicated, not the simple linear depen- 
dence one might have imagined. Instead, the dependence 
on the azimuthal angle is different at a different radial 
coordinate, indicating a violation of the adiabatic condi- 
tion with respect to the radial motional state. However, 
as long as the internal state can follow adiabatically the 
ffipping z-bias B-field, our intuitive vortex pump proto- 
col remains effective. As we demonstrate in Fig. [3fa), for 
successive ffipping times of 30 ms to 60, 90, and 120 ms, 
the phase distribution approaches the same form with 
the increasing level of adiabaticity. 

Next, we study the effect of misalignment of the ax- 
ial symmetric optical trap with respect to the geometric 
center of the external B-field. Any misalignment will 
reduce the axial symmetry, and thus break the angular 
momentum conservation law [40] and adversely impact 
the quality of our vortex pump protocol. Without loss of 
generality, we assume the optical plug's center is (xq, 0, 0) 
while other symmetries are assumed to remain. After 
extensive numerical simulations as before, we find that a 
misalignment of the percent level of the transverse size 
for the condensate is already detrimental to the pump 
protocol as shown in the upper row of Fig. ^h) when 
the bias fiip time is 30 ms. More generally, we find that 
the larger is xq, the stronger is the deviation from the 
intended vortical phase distribution pattern. Somewhat 
surprisingly, however, we find a simple solution to this 
ailment. By increasing the z-bias fiip time, the compli- 
cated phase structures disappear gradually as shown in 
the lower row of Fig. [3l[b), where the fiip time is taken 
to be 60, 90, 120 and 150 ms from left to right. These 
results show that unless there is a total failure in aligning 
the optical trap with respect to the B-field trap, the nom- 
inally small misalignment can be compensated for by the 
proposed vortex pump protocol with a correspondingly 
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FIG. 3: (Color online). The condensate phase distribution for 
the Mf — 1 state after flipping the 2;-bias in an IPT. (a) The 
ideal case of an IPT with perfect alignment with the optical 
trap and the plug. The flip times are 30, 60, 90, and 120 ms 
from the left to the right, accompanied by an increasing level 
of adiabaticity. (b) The same as in (a) except for a misaligned 
optical plug Vp^ = [/exp(-((x - xof + y'^)/2pl), U/Ti = 
2 X 10^ Hz, and po = 5 jim. From the left to the right in the 
upper row the results are compared for xo/po = 0.05, 0.1, 0.2, 
and 0.3, respectively, at a bias flip time of 30 ms. In the lower 
row xq — 0.2y9, where results from different flip times of 60, 90, 
120, and 150 ms are compared from the left to the right, (c) 
The same as in (a) except for a non-axial-symmetric optical 
trap Vo^ Mu\[(l + £)^x^ + (1 - eYy'^]/2 parameterized 
by e. The upper row is for e — 0.05, and e = 0.1 in the lower 
row. The z-bias flipping times used are 60, 90, 120, and 210 
ms from left to right, respectively. 



slower fiipping time. 

Finally, we study the potential degradation due to an 
asymmetric optical trap. As with the misalignment con- 
sidered above, an asymmetric optical trap will void the 
conservation of Dz . Although the intended adiabatic ma- 
nipulation of the B-field can still transform all atoms from 
the Mp = — 1 to Mp = 1 spin state and cause the Fz to 
change from — to for a spin-1 condensate, the axial 
motional angular momentum per atom gained will not 
be equal to 2fi, except when the asymmetry is small. In- 
stead, Lz generally is found to oscillate near the end of 
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FIG. 4: (Color online). The motional angular momentum 
gained with our protocol in an asymmetric optical trap and 
an IPT. (a) The temporal development of Lz during the flip 
of the z-bias field at an asymmetry parameter e — 0.05. The 
flip time T is 30, 60, 90, 120, and 150 ms, respectively, and is 
shown by off-set lines from the lower to the higher. The off-set 
shift for each curve is given by — 0.25/ix [(T/Tm)-1], and 
Tm — 30 ms. (b) The difference between angular momentum 
gained from our protocol in an asymmetric optical trap with 
= 2% predicted for the ideal case. The curves denote 
different asymmetry parameter e from 0.05 to 0.15 as shown 
by the insert labels, (c) The same as in (b), but now for the 
e dependence at the same flip time T. The five curves are 
labelled by their respective T /Tm values from 1 to 5. 



the B-field flipping, as we show in Fig. [H For relatively 
small trap anisotropy, for instance at £ = 0.05, the vortex 
generated is not far away from the ideal case of perfect 
symmetry. This can be readily derived from the axial an- 
gular momentum shown in Fig. [4](a) or by comparing the 
condensate's phase structure in Fig. [31(c) with Fig. [31(a). 
With larger optical trap asymmetries, for example at the 
threshold of £ = 0.1 for the system and parameters that 
we consider, the average angular momentum per parti- 
cle {Lz) gained by the condensate is no longer an integer 



multiple of Ti. The condensate phase structure becomes 
more complicated at short flip times, and if the optical 
trap asymmetry is distorted further a dynamical insta- 
bility begins to set in within the same evolution time. 
Based on our numerical simulations, when e = 0.125 and 
the flipping time is set to 30 ms, a dynamical instability 
occurs. In addition, we find significant blurring in the 
condensate density distribution after flipping the B-field. 
The phase difference between nearby spatial positions is 
totally smeared out. As the flipping time is increased 
further, dynamical instability or the blurring disappears 
while high quality winding pattern reemerges. Our simu- 
lations also show that the dynamical instability survives 
at larger rotations for smaller trap asymmetries. 

III. GEOMETRICAL PROPERTIES OF STATIC 
MAGNETIC FIELDS 

According to Maxwell's equation, the spatial distribu- 
tion of a static B-field can be expressed as the gradient 
of a scalar potential ^(r), i.e., 

B{f) = [B,{r), By{r), B,{r)] = V^- (8) 

Since the divergence of a static B-field is always zero, 
or V • 5 = 0, the scalar potential satisfies the laplace 
equation V^<l>(r) = 0. Thus $(r) can be expanded in 
terms of the spherical harmonic functions, to 

^(^) = ^cimr^Pf^ {cos 0)e''^'^, (9) 

Im 

in the spherical coordinate (r, 6>, (j)). 

The vortex pump protocol we proposed relies on the 
repeated flips of the z-bias field adiabatically from 
the -\-z to the —z direction and vice versa. The vorticity 
gained by a condensate can be understood in terms of 
the conservations of Dz or Lz -\- Fz, which are ultimately 
determined by the angle between B{r) and the x-axis 
(or the 7/- axis). A vortex state has a two-dimensional 
phase structure. For a stable vortex state with its angular 
momentum pointed along the z-axis, the velocity v = 
^V(j)p must flow along the same azimuthal direction, and 
the phase in any plane with z = Zq must be similar 
to that in the z = plane, i.e., z^ (j)) = 4'p{p^ 0, (j)) + 
f{z). To effectively obtain a stable vortex state, the angle 
1^ should (i) depend only on (j) but remain independent 
of i.e., ^{p^z^cj)) = 0,^); or (ii) depend on both 
(j) and z but take the form i^{p, (j)) = i9{p^ 0, (j)) + g{z). 
These constrains then lead to the condition 

Clm = 0, (10) 

where m 7^ ±/ or m 7^ on the expansion coefficients. 
To obtain a vortex state with a definite winding number, 
the scalar potential thus has to take one of the following 
two possible forms 

^[^\f) (X pUos{l(j)^r]), (11) 



or 



(12) 



where (p, ^) are the cyUndrical coordinates and 77 is a 
constant angle. In the derivation of the second form, we 
have used the fact that ^(f) is real. 

In the first case above and as was discussed in the thesis 
work of A. E. Leanhardt [39], the 5-field takes the form 



cos{l — l)(j)ex — sm.il — 1) 



(13) 



with I > 2. The corresponding conserved quantity be- 
comes Dz = Lz — {I — I) X Fz. When / = 2, this is the 
field of an IPT. In the work of Mikko Mottonen [38] , both 



fields B^\f) of a hexapole and B2^{r) of a quadrupole 
are used. For a spin-1 condensate, the vorticity gained is 
Lz = 2{l — l)h after adiabatically flipping the Bz{r). 
In the second case, the simplest B-field is 



B^''\r)=B,{r)e,^Bz{r)ez 



(14) 



as in a 3D QT. In the x-y plane, B^^^"^ points along the 
direction of ip. For a spin-1 condensate, the vorticity 
gained is Lz = —2?i per atom after the flipping of Bz{r). 
For completeness, we have confirmed this numerically for 
a 3D quadrupole trap. A slight complication arises in 
this case, as the B-field depends on the z-coordinate, 
and its zero value point moves from —00 to +00 dur- 
ing the flipping of the z-bias from —0.1 to 0.1 Gauss. 
To demonstrate the numerics, we used the coupled 3D 
Gross-Pitaevskii equations. The gain of a = —2h vor- 
tex also can be understood in this case as being due to 
the conservation of Fz + Lz [4Q] . 



IV. THE OPTIMAL VORTEX PUMP 
PROTOCOL 

As discussed above in some detail, starting in the ini- 
tial spin state |(— 1));^, vorticity can be pumped into the 
motional state of a condensate if the z-bias field Bz is 
flipped from the positive to the negative z-axis. When 
the static field is of the form B\^\f)^ the vortex state 
has L z = 2(1 — 1)%] li the form is changed to B^^^\r)^ a 
vortex state of Lz = —2fi is obtained. The geometrical 
constraint on a static B-field together with its coupling to 
the atomic hyperfine spin stops the change of the atom's 
internal state through the flipping of the z-bias B-field. 
In fact the opposite vortex state oi—2(l — l)?i with / > 3 
can never be created through the flipping of Bz associated 
with any B-field. Clearly, our studies provide an optimal 
vorticity pumping protocol based on a cyclic flipping of 
the bias Bz field as devised below. 

If we have access only to the field b\^^ (r) , we can turn 
it on and generate a vortex with Lz = 2{l — l)h through 
the flipping of the bias in the first step. Then, this is 




FIG. 5: (Color online). Different vortex pump protocols com- 
pared: from the top to the bottom. The parallelogram and 
hexagon denote, respectively, 2D quadrupole and hexapole 
B-fields, and the stacked pyramids denote a 3D quadrupole 
B-field. The initial z-bias is denoted by the solid black line 
arrow that is flipped to the opposite direction denoted by gray 
line arrow in the first step, the slanted solid black line arrow 
in the right figure of (c) refers to the transverse bias, (a) The 
original theoretical protocol of Ref. [2511 , which has been ex- 
perimentally demonstrated [s^, 1351 , l36ll can pump only once, 
(b) The continuous pump protocol of Ref. [3^]. It works, but 
the condensate vorticity does not monotonically change; (c) 
and (d) are the optimal protocols we discuss. They are capa- 
ble of continuous and monotonic changes to the vorticity of 
a condensate. They are the optimal protocols based on the 
manipulation of B-fields. 



followed by a second step where the field B^^^ (r) is turned 
off, and the bias field Bz is flipped back with the help of a 
constant transverse bias B^x as we have suggested. The 
second step does not create any vorticity but returns the 
internal state to what it was initially. Repeating the two 
step protocols, we pump Lz = 2{l — l)h units of vorticity 
into the condensate for each cycle. The motional angular 
momentum gained from this protocol for the first five 
repetitions is shown in Fig. [6^a). 

If we have access to both fields B\^\f) and B^^^\r)^ 
the second step of each cycle is modified by tuning on the 
field B^^^\r), which gains an additional 2?i units of vor- 
ticity in each cycle, for a total of 21% angular momentum 
gained per cycle. Figure [61(b) shows the continuously in- 
creasing vorticity for this more general case of the two 
type 2D multipole fields with / = 3 and 1 = 2. 

Finally, in light of our analysis above, the protocol de- 
vised in Ref. [s^, which uses only the field of B[^\f) 
with / = 3 and / = 2, respectively, in the first and the 
second steps of each cycle, is clearly not optimal. In the 
first step, the condensate gains a vorticity of 4fi, while 
in the second step, it gains a vorticity of —2fi^ or loses 
2fi. In the end, the net gain is only a 2fi per cycle. We 
summarize the various vortex pump protocols based on 
the manipulations with external B-fields in Fig. [5l 

In actual experiments, many complications could arise 
that make the selection of the optimal vortex pump pro- 
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tocol an entirely system-dependent matter. For instance, 
the loss of atoms and quantum coherence due to dissi- 
pation or decoherence strongly limit the lifetimes of all 
currently available atomic superfluids. If the criterion for 
the optimal protocol is defined as the maximal amount of 
vorticity gained during a fixed time, our analysis above 
then points to a clear winner: the repeated applications 
with only the hexapole field ^3^^ (r) as evidenced by the 
results in the illustrated figures above. 

Before concluding, we want to point out that Dz re- 
mains conserved within our formulation using the ex- 
act eigen-energies for the three Zeeman states, which 
is indeed confirmed in numerical simulations. This is 
easily understood because the effective Zeeman inter- 
action H^^ = Eo - r]oF ■ B/B + 5{F ■ Bf/B^ with 
7/0 = (£^-1 - Ei)/2 and 8 = {Ei ^ E_i - 2Eq)/2 com- 
mutes with Dz. In addition, Lz is also conserved during 
the 10 ms intervals when only a uniform B-field is present. 



V. CONCLUSION 

In conclusion, we have proposed and demonstrated nu- 
merically a simple yet efficient vortex pump protocol ca- 
pable of continuously increasing or decreasing the vor- 
ticity of a condensate through repeated manipulations of 
external B-fields: a bias along the axial direction, a 2D 
quadrupole field or hexapole field, and a 3D quadrupole 
field or a transverse bias field. Based on a general consid- 
eration of the geometrical properties of static B-fields, we 
have shown that the above choices of B-fields are the only 
possibilities capable of continuously controlling angular 
momentum vorticity in a spinor condensate through the 
manipulation of B-fields. In Fig. [6] we compare the two 
optimal protocols for continuously increasing vorticity for 
the first five repetitions. The vorticity is seen increas- 
ing to 20% and lOfi by using only B[^\r) for a hexapole 
(/ = 3) and a quadrupole (/ = 2), respectively, in (a), and 

to 30fi and 20fi by combining both b[^^ (r) and B^^^^ (r) 
for / = 3 and / = 2, respectively, in (b). To assure adia- 
baticity throughout the B-field manipulations, each cycle 
has to be engineered differently in the two protocols. For 
the first protocol as shown in (a), the duration for each 
cycle, or its period is Ti = 40 ms. During the first 30 ms 
we use B[^\r) to increase vorticity, and in the following 
10 ms the transverse bias along the x-axis is turned on to 
reset to the initial internal state through a second z-bias 
fiip. For the second protocol of (b), the period required 
becomes T2 = 60 ms. During the first 30 ms B[^\r) is 
present and in the following 30 ms B^^^^ (r) is turned on 
while b[^^ (r) is turned off. 

Our discussions indicate that for / = 2, the second 
protocol is more efficient if a maximal amount of vorticity 
is sought. While for / = 3, our results show that the 
efficiencies for the two protocols are about the same due 
to the different time requirements for adiabaticity. Based 
on the experience from the numerical simulations and 
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FIG. 6: (Color online). The continuously increasing vortic- 
ity of a condensate for the optimal protocols, (a) with only 
Bi^\f) for a hexapole (/ = 3 the upper blue curve) and a 
quadrupole (/ = 2 the lower red curve). Ti = 40 ms is the 
time of flipping. In the first 30 ms Bi^"^ (r) is turned on to gain 
vorticity, while no vorticity is gained in the following 10 ms 
when a transverse bias along the x-axis is turned on to reset 
the internal states to be as they were initially, (b) with both 
B-fields Bi^^ (r) (/ = 3 for the upper blue curve and / = 2 for 
the lower red curve) and B^'^^^f) available and used for the 
first and second steps, respectively, of each cycle. T2 = 60 
ms is the total time for each cycle. In the first 30 ms B^^"^ (f) 
is present and in the next 30 ms B^^^\r) is turned on while 
B[^\r) is turned off. The vorticity increases at all times in 
this case. The dashed lines in two figures are guides for the 
eye. 



from experimental works, it seems more practical to use 
the first protocol. The reason is rather simple, it is always 
more difficult to have two types of magnetic traps (from 
two sets of current carrying coils) working together. The 
first protocol requires only one type of multipole 2D B- 
field B[^^\r), that can be realized on an atom ch ip |43 |, 
using a three- or four- wire structure for a 2D QT [45] , or 
a five-wire structure for a 2D hexapole trap [46|] . 
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We demonstrate above the practical limit of reaching 
30h per atom using our protocol, which is already much 
larger than anything reported experimentally through 
the stirring of a condensate. We do not know what the 
maximum limit is. In fact, the exact value of this maxi- 
mum is perhaps not so important either for two reasons: 
1), our theory is mean field and may well break down at 
such larger angular momentum per atom because quan- 
tum correlations become essential; 2), we are also limited 
by the computational resources available for such large 
scale simulations. 

The prospect for realizing the proposed vortex pump 
mechanism provides new impetus to active pursuits of 
quantum simulations of strongly interacting_many body 
electronic systems in terms of cold atoms [3, 4, [H, 0, 
i, i, mi, m, [H, Q m, EB] . it opens new avenues of 
theoretical and experimental research into the coupling 
between the internal and motional degrees of freedom and 
shines new light on relevant topics in spintronics studies 

Finally we note that in all our protocols discussed 
above, there always exist a conserved quantity Dz for the 
ideal case of the model system. To create a vortex, one 
simply needs to adiabatically fiip the B-field to transfer 
all atoms from the weak field seeking Zeeman state | — 1)^ 
to I + 1)2, or the reverse. An obvious connection can be 
made with the famous Einstein-de Haas effect to deepen 
the understanding of our protocol and to possibly extend 
to the systems of dipolar quantum gases [48*, . During 
our extensive numerical simulations, we find that a small 
B-field gradient or 5^, combined with a stronger opti- 
cal plug enforces the stability of the high vorticity states 



as generated. As is well known high-order multipole B- 
field traps are usually not ideal candidates for magnetic 
traps because they are too weak to trap atoms. As far 
as vortex pump is concerned, however, a hexapole trap is 
preferred over a quadrupole trap because we are making 
use the topological structure of the B-field, rather than 
its spatial dependence when trapping is considered. 

To help with further experimental effort, we suggest a 
feasible implementation scheme on an atom chip. As is 
well known, an IPT can be realized by combing a chip 
based Z-type trap with a constant B-field [i^] • This type 
of IPT is not an ideal setup for implementing our proto- 
cols because the locations of trap zeros at the vanishing 
B-field change following increasing or decreasing currents 
in the on chip wires. The IPTs from the more stable 
three- or four-wire structures, on the other hand, are 
more promising candidates as evidenced by their demon- 
strated successes and stabilities in our chip matter wave 
manipulations [45.]. For the three- wire structure, three 
equally spaced parallel wires lay on the surface of an 
atom chip, with the current in the middle wire oppo- 
site to the other two side wires. The four- wire struc- 
tures are analogously configured. Following similar ideas, 
a two-dimensional hexapole trap can be constructed by 
five wires [46|. Additionally these chip-based traps can 
be combined with optical traps, such as an optical plug, 
to realize a rich variety of confinement geometries. The 
force of gravity can be made to align along any of the 
symmetry axes of the two-dimensional multipole trap or 
compensated with an additional optical confinement. 

This work is supported by US NSF, and NSFC and 
MOST of China. C.R. acknowledges support from ARO. 
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